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ABSTRACT 

Production of large-pr hadron pairs by polarized photon on a longitudinally polarized pro- 
ton towards probing the polarized gluon distribution is studied. Resolved photon contributions 
and the effect of changing the scales are taken into account, and predictions are presented. 
A very recent experimental result at cm. energy 7.18 GeV is compared to our predictions 
extented down to this energy. A proper combination of cross sections is also considered. 
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I. INTRODUCTION 



The determination of the size and shape of the polarized gluon distribution Af? remains 
a major problem in Spin Physics. Clearly, the way to proceed is to study theoretically and 
experimentally polarized reactions dominated by subprocesses with gluons in the initial state. 
To this effect, experiments on charm production by polarized photons on longitudinally po- 
larized protons (polarized photoproduction) [1], large-pr direct photon and jet production in 
polarized p-p collisions [2], etc. will be carried [3,4]. 

As a reaction leading to useful information, Ref. [5] has proposed the production of large- 
Pt hadron pairs Hi, H2 in polarized photoproduction: 

^+p- ^ H1 + H2 + X (1.1) 

An experiment could well be carried in COMPASS. 

In view of this, we have undertaken an independent study of the reaction (1.1). Working, 
as in [5], at Born level (leading order (LO) in a^) we differ in the following from [5]: 

(i) We take into account the resolved photon contributions, which are left out in [5] . 

(ii) In general, reaction (1.1) is dominated by the subprocesses 

(a) — > qq, (b) 7" — > qg (1.2) 

In [5], the first is well taken into account, but the second is treated in a rather unclear 
way. Here the subprocess (b) is treated on equal footing with (a). 

(iii) We consider the effect of changing the renormalization and factorization scales; in [5] 
this effect has also been left out. 

(iv) In [5] the fragmentation of the final partons to hadrons is treated via Monte-Carlo 
methods, which somewhat obscure the procedure. Here we use the conventional QCD 
approach with recent fragmentation functions [6]. 

(v) Very recently Ref. [7] presented an experimental result on (1.1); this is discussed and 
compared to our predictions. 

(vi) We show that a proper combination of cross sections for certain choices of Hi and H2 
will make a more clean probe. The combination, however, involves four cross sections, 
and the experiment will be more difficult. 

Furthermore, apart from the cross section calculated in [5] and in relation with [7] (Adcr/d^idx, 
for the definition of 0i and x see Sect. II), we present also results for the transverse momentum 
distribution Ada/d4>idxT- 



With the COMPASS experiment in mind (polarized muon-proton scattering), we take into 
account that the (initial) photons are in general quasi-real (7*). 

Sect. II presents our general formahsm for the cross section Ada/d(j)idx and Sect. Ill for 
Ada/d(f)idxT- Sect. IV presents results for Ada/d(f)idx and the corresponding asymmetries. 
Sect. V presents results for Ada / d(f)idxT- Sect. VI presents the above mentioned combination 
of cross sections as well as our results. Finally, Sect. VII presents our concluding remarks. 



II. GENERAL FORMALISM FOR MG/d(j)idx 



The reaction (1.1) has, to some extent, been studied in Ref. [8], and here we avoid 
repetition as much as possible. Consider the contribution of the subprocess 

t (pi) + t (P2) Ci (ps) + C2, (2.1) 
where the quantities in parentheses denote 4-momenta, and let 

S=(pi+P2)% t={p^-pif, U={p3-P2f (2.2) 

{s + t + u — 0). Neglecting intrinsic transverse momenta, the hadrons Hi, {i — 1,2) are 
produced in opposite hemispheres with transverse momenta kiT and cm. pseudorapidities rji 
with respect to the photon. Denoting by y/S the total cm. energy and by 0i the azimuthal 
angle of Hi and introducing 

XiT = 2kiT/ VS, 

it follows that the cross section for (1.1) is formally given by [8,9] 

d<PidxiXldrjidrj, = f / ^"^^^"/^ ^""^ ' ^'-'^ 



where AFf,/^ the polarized momentum distribution of parton b inside the photon and 

Aa = ^AFa/p {xa) ^%Dh,/c, (zi) Dh,/c, {Z2) ; (2.4) 

the limits of integration in (2.3) are specified later. In (2.4), Ada/dt is the cross section for 
the subprocess (2.1), DHi/a (^i) is the fragmentation function for Cj H^ and 

Xa = Xb exp {-r]i - ri2) (2.5) 



Zi = XiT (exp (?7i) + exp (772)) f^Xb 



(2.6) 



Eq. (2.4) expresses the contribution to the physical cross section from both direct and resolved 
7, the former corresponding to AF^/^ {x) = S {1 — x). The cross sections for the subprocesses 
(1.2) are: 

~W ~ s2 ~W ~ 3s2 -St ^ ' ' 

The corresponding cross sections for the resolved 7 contributions are taken from [10] with 
t ^ u (see also [8]). 

Here we define the variable 

X = exp (-771 - 772) (2.8) 

and determine first H^da / d(f)idx . Also introduce 

h = exp (772) (2.9) 

Taking, as in [8] , the x-z plane to be defined by p2 and ki (i.e. 0i = 0) we may write still in 
a formal way 

(5", 01 = Q^x) = J dXiT j dX2T J dXbAFh/^ (Xb) J "^^Cr {S,X,Xb,XiT:X2T,h) 



#1 



(2.10) 



The physical meaning of the variable x is clear from Eq. (2.5): it for Xb — 1 (direct 

7). 

The limits on h are specified by the condition Zi < 1, which, in view of (2.6), (2.8) and 
(2.9), implies: 

h + x~^hr^ < Xxb, 
where A = min {2/xit, '^/x2t)- We find: 

h- <h < h+, 

where 

4\ 1/2" 



Xxb ± (^X^xl - -) 



h-h+^l/x (2.11) 



Clearly, we must have 

Xxb > 2/v^ (2.12) 
Denoting the lower limit of Xit, X2t integrations by x^"* (= 2ki^^ /y/S, k^^ to be fixed by 



experiment) we write 

/ dX2T = / dX2T + / dX2T (2.13) 



So we find: 



dcpidx 4 Jx%^ 1^ 



where 

f-Xbv^ fXxT [■h+ 



rxbyjx rxiT rri+ (III 

/i = /^^ dxir L^^ dx2T —AFb/^{xb)Aa{S,x,Xb,xiT,X2T,h) 

cc-j-T Xrp J h It/ 



with 
and 

f^bV^ rx2T rh+ (ill 

h = I dx2T I <ia;iT / —AFb/^ {xb) Aa {S, x, Xb, xit, X2t, h) 

with 

/i-j- ^— X2^Xb i ( ^2T^b ^ / ^ 1 

Given a;^'*, the condition (2.12) determines the minimum value of x allowable. 

X > {xj} /Xb^ 

and since Xb <V. 

M 



III. GENERAL FORMALISM FOR Ada/d(/)idxT 

We start again from Eqs. (2.3), (2.4) and change variables: 

1 

X2T ^ XT = - {XIT + X2t) , Tji ^ hi = 1 = 1,2 



so 

Ada 



d(f)\dxTdxiTdh\dh2 2h 
and Act given by (2.4). Now 

Xa = Xb/hih2 



S f 

— / dxbAFb/j{xb) Aa{S,XT,Xb,XiT,hi,h2) 

lAl2 J 



and 

^1 = ^ {hi + h2) , Z2 = — [hi + h2) 

The conditions Zi < 1 and < 1 imply 

where A = min {2/xit, 2/ {2xt — Xit))- As in Sect. II: 

where now 



1 



Axfo ± {}?x1 — Axb 



1/2 



h-h+ = Xb 



Here we have the condition 



Xb > 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



Clearly, for xit < xt' X — 2/ {2xt — a^ir), whereas for xyr > xt' A = 2/xit- So in the 
present case we write 



The final result is 



/ dXiT = / dXiT + / dXiT 



(3.9) 



(3.10) 



where 



Ji = / dx,T / dxb ^ (^f) (3-11) 

Jx^j, Jh- n2 Jxb/h2 ill 

with X2T = 2xt — XiT and 

XlT,min = max (a;5?\ 2xt - l) , h^rnax = '2Xb/x2T - /i2 (3.12) 



h 



^^^±{^-Xb 
X2T \^2T J 



1/2 



and 



rXlT,max 1-1 rh+ dfl^ rhl^max dkl . ^ / x a 

J2 = / dxiT dXb — —AFb/^ {Xb) Act 

Jxt ''^2T ""2 Jxb/h-i 111 



(3.13) 



(3.14) 



with 



2^ir,max = min (2xt - \ l) , hi^rnax ^ "^Xb/xir - h2 (3.15) 

= — ± f ^ - x.y^' (3.16) 

In determining XiT^min wc took into account that Xit — 2xt—X2t > 2xt— 1, and in determining 
xiT,max that Xit = 2xr - X2t < 2xr - xP . 

IV. RESULTS FOR Ma / d(j)idx AND THE CORRESPONDING 

ASYMMETRIES 

We present results for the three sets A,B,C of LO polarized distributions of [11], which 
can be roughly characterized as follows in terms of A^f [x] (= ^Fg/p [x, Qo)): 

Set A: Agf (x) > and relatively large 

Set B: Agr {x) > and small 

Set C: Agi (x) changing sign; Ag( [x) < for a; > 0.1. 

The fragmentation functions -D^./cj are taken from [6] (LO sets). In (Q) we use A = 0.2 
GeV and 4 flavors. The rcnormalization and factorization scales arc taken equal and with a 
central value Q = Qc = kir + k2T- We first present results at a typical COMPASS energy 
= = 12 GeV and for 4°^ = 1.4 GeV [5]. 

Regarding the resolved 7 contributions, we have used the maximal and minimal saturation 
sets of the polarized photon distribution functions of [12,13]. We have also carried calculations 
with the distribution functions of [14], belonging to the class of the so-called asymptotic 
solutions, and we simply report the results. 

To account for the fact that the photons are quasi-real we multiply (2.4) by the Weiszaecker- 
WiUiams factor: 

A/ iy) = —AP^i (y) /n ^^°^ (4.1) 

where AP^^ (y) = 1 ^ (1 ^ v)^ IVi = muon mass, and we take a typical value Qln^x ~ ^ 
GeV^\ for incident lepton cm. energy (corresponding to Ei = 200 GeV [5]): y = S/Si. 

Fig. la presents Ada/d(f)idx for direct and resolved 7* contributions with = tt^ or 7r~ 
(fragmentation functions (A.4)-(A.8) of [6]). The presented resolved contributions correspond 
to the maximal saturation set of [12,13]; those of the minimal are somewhat smaller. So, in 
general, the resolved contributions are much smaller than the direct. However, in particular 



for the set A of [11] and in the range 0.15 < x < 0.2, where the direct contributions change 
sign, the resolved are not insignificant. The asymptotic solution of [14] gives even larger 
resolved contributions. 

Notice that the differential cross sections for the direct 7 contributions change sign at 
some X < 0.2; this is due to the two competing subprocesses of Eq. (1.2). At the lower x, (a) 
dominates, whereas at higher x, (b) takes over. Hence the place to obtain information about 
AG is at the lower first was pointed out in [5]. 

Fig. lb presents the asymmetries 

^ ^ Ada/d(j)idx 
da/d(f)idx 

for the sum direct+resolved and again Hi — tt"^ or tt" . For the unpolarized da/d(l)idx we use 
the CTEQ distributions [15] and the photon distribution functions of [16], LO sets. Here, to 
account for quasi-real photons, in (4.1) we replace AP^i by P^i (y) — 1 + (1 — y)^ /y; hence 
the asymetry is reduced. For each of the sets A, B and C the strong line corresponds to the 
central value Qc = ki'r + k2T- For the sets A and C, Fig. lb presents also the effect of changing 
the scales in the range Qc/^ < Q < 2Qc. 

Fig. lb also presents an estimate of experimental errors using the expression 

= ] (4.3) 

We take beam polarization Pb — 80%, target polarization Pt — 25%, pion-kaon detection 
efficiency e — 1 and integrated luminosity L — 2 fb~^ [5]; in (4.3) a^*p is the unpolarized cross 
section for quasi-real photon- proton scattering integrated over a bin Ax — 0.17. 

On the basis of Fig. lb we conclude the following on the experiment: First, the sets A 
and B cannot be distinguished. Second, the sets A and C can barely be distinguished in the 
small range 0.15 < x < 0.2. We note that at smaller x the cross sections 7*p become much 
smaller and 6A^*p much larger. 

Now we turn to kaon production, and Fig. 2a presents Ada/dcpidx for = or K~ 
(fragmentation functions (A.19)-(A.23) of [6]) and Q = Qc- The presented resolved contribu- 
tions are as in Fig. la; now for both sets A and B of [11], at x > 0.25, they are important. 
Fig. 2b presents the corresponding asymmetries together with the effect of changing the scales 
and an estimate of the experimental errors, as for Fig. lb. Now the latter are significantly 
larger (smaller cross sections), making very difficult the distinction even between sets A and 
C. 

In [5], apart from kaons, the production of charged hadron pairs is considered. The unpo- 
larized cross sections for the production of charged hadrons are, of course, greater than those 
of charged pions only. Thus the estimated errors will be somewhat smaller. 

As it has been stated, very recently an experimental result was presented for (1.1). Its 
energy is low, \/Si — 7.18 GeV , and therefore kiT limited; also, the way one reaches the final 



result is somewhat unclear. Nevertheless, in view of its importance and of the fact that it is 
the first experimental result, it is perhaps of interest to extend the calculation of Sects II and 
IV down to y/Sl = 7.18 GeV. This experiment selects events containing at least one positively 
charged hadron and at least one negatively charged hadron. Hence the fragmentation functions 
sould be separated to those for 7r+ {K'^) and for tt" {K~). For the separation we use the 
subsequently presented expressions (6.3)-(6.6). The fact that the experiment was carried at a 
low energy necessitates the choice A;^ -* = 1.1 GeV. Taking into account also the virtual photon 
depolarization factor D = 0.93 [7], the predicted asymmetries together with the experimental 
result are shown in Fig 3; clearly, Set A (or B) is favored. 

It is interesting also to note that the effect on the asymmetry of changing the scales is 
small. 

The values \^ — 12 GeV and k^^ = 1.4 GeV imply certain limits on the rapidities rji and 
the invariant mass of the hadron pairs m{HiH2), which amount to acceptance cuts. Taking 
the variable x in the range 0.055 < x < 0.8, the integration limits for the variable h in 
Eqs (2.15) and (2.17) combined with Eqs. (2.9) and (2.8) imply the fohowing hmits on the 
rapidities: -1.9 < % < 2.0 and -1.7 < 772 < 2.1. These limits imply m{HiH2) > 2.81 GeV, 
which is in accord with [5]. 

It should be remarked that, instead of = kiT + k2T, the choice Qc = {kiT + k2T) /2 is 
also reasonable. Then varying the scales in the range (5c/2 < Q < 2(5c, near the lower limit, 
with kj^ — 1.4 GeV, we enter a region where perturbative QCD is uncontrollable. One cannot 
take kj^ much larger because Ada / d(l)idx becomes too small. 

V. RESULTS FOR Ada/d(l)idxT 

The transverse momentum distributions Ada / d(f)idxT and da/d(f)idxT are calculated for 
the same distributions and fragmentation functions as Sect. IV, as well as for \fS — 12 
GeV and k'^^ = 1.4 GeV. We present results only for Q — Qc = kiT + k2T as functions of 

XT = (kiT + ^2t) / V^- 

The indicated errors have been estimated on the basis of Eq. (4.3) with the unpolarized 
cross section integrated over a bin in corresponding to ApT = 1 GeV. 

Fig. 4a presents asymmetries for Hi = 7r+ or tt~ . Clearly, even without accounting for the 
variation of the scales, sets A and C, as well, are hard to distinguish. 

Fig. 4b presents asymmetries for Hi — K'^ or K~ . The conclusions are the same as for 
Fig. 4a. 

Again, as in Sect. IV, taking the variable xt in the range 0.25 < xt < 0.8 we obtain the 
following limits: -1.9 < 771 < 2.1, -2.1 < 772 < 2.1 and m{HiH2) > 2.79 GeV. 



VI. THE COMBINATIONS OF CROSS SECTIONS 



Denote, for simplicity, a {H1H2) either of the cross sections Ada / d(f)idx and Ada / d4>idxT 
for ^ + p*' — )• Hi + H2 + X. As it is discussed in Refs. [17] and [8], neglecting the resolved 
7 contributions, the combinations 

A{tt) — a {t:^t:~^ + a {t:~t:^^ — a {tt^t:^^ — a {t:~t:~^ (6.1) 

and 

A{K)^a (k+K-) + a (k-K+) - a (k+K+) - a (r- R-) (6.2) 

isolate the contribution of the subprocess — > qq. 

When the resolved 7 contributions, calculated via the polarized distribution functions of 
[12,13] or [14] are taken into account, the contribution of — > qg is not completely 
eliminated, but we find that for x < 0.4 it is smaller by ~ 2 orders of magnitude than the 
contribution of 7^^ — > qq. Hence the difference in Ag between the sets A, B and C is 
displayed much better. Below we present results for the corresponding asymmetries and for 
Q = Qc = kiT + k2T. 

The calculation of (6.1) requires the separation of the fragmentation functions for tt'^ and 
7r~. To this purpose, as in [8], we use: 

I + z 

Dn+/u {z) /D^-ju {z) = D^-id {z) jD^+ia {z) = —— (6.3) 

-L Z 

and 

{z) /D,-,g {z) = D,+,s {z) ID,- 1, {z) = 1 (6.4) 
For the calculation of (6.2) we use [8]: 

1 + 2; 

Dk^Iu {z) /Dk-/u {z) = Dk+/s {z) /Dk-/s {z) = —— (6.5) 

-L Z 

and 

DK+/g {z) /Dk^/s {z) = Dx+fd {z) /Dx-fd (z) = 1 (6.6) 

The effect of changing the scales is very similar to that of Figs lb and 2b. 

The indicated errors have been estimated as follows: First, for each of the cross sections 
a {H1H2) we have determined an error (via Eq. (4.3)) and then we have taken the square root 
of the sum of the squares of these errors (assuming independent measurments of a {H1H2)). 

Fig. 5 displays asymmetries corresponding to cross sections Ada / d(j)idx] Fig. 5a refers to 
A (tt) and 5b to A [K] . Now the differences between the sets A, B and in particular C are 
larger and over a wider range of x than in Figs, lb and 2b. The errors are larger, but A (tt) 



appears to be useful in distinguishing between sets A and C. As for A (K), the errors are too 
large to be of any use. 

Fig. 6 displays asymmetries corresponding to cross sections Ada / dcpidxT] Fig. 6a refers to 
A (tt) and 6b to A (K). Again, the differences between sets A, B and in particular C are larger. 
Fig. 6a seems to show that A (tt) does distinguish between sets A and C at 0.25 < x < 0.3. 
Again, for A (K) the errors are too large. 

Of course, as in Sects IV and V, we present predictions for tt and K separetely. In an 
experiment detecting tt + K the errors will be somewhat smaller. 



VII. CONCLUDING REMARKS 



On the basis of Fig lb we have concluded that the sets A and C of polarized parton 
disrtibutions can barely be distinguished and only in the small range 0.15 < x < 0.2. The 
sets A and B cannot be distinguished. 

Nevertheless, as we stated, in an experiment detecting all charged particles, the errors will 
be smaller. Furthermore, one way to increase the asymmetry (4.2) is by increasing — ^J~Syp] 
this will somewhat reduce the implications of producing quasi-real instead of real photons. 
Alternatively, one may decrease the incident lepton cm. energy ^JSip- On the other hand, pro- 
ducing quasi-real photons with ^fS very near ^JSip might make the experiment more difficult 
[8]. 

Polarized real photons at energy comparable to that of the present paper \fS (~ 10 GeV) 
are available at SLAG. 

A very recent experiment at \fS = 7.18 GeV [7] favors set A or B (Fig. 3). However, 
the fact that this energy is low and the way the final result is obtained makes necessary the 
repetition of the experiment at a higher energy [1] as well as experiments at even higher 
energies with different reactions involving polarized initial particles [2]. 

A somewhat better probe of Agr appears at first sight to be the combination A (tt) of 
cross sections corresponding to Ada/d(l)idx (Fig. 5a) and even better the combination A (tt) 
corresponding to Ada/d(f)idxT (Fig. 6a). However, in our estimate of errors, only statistical 
ones are taken into account. The systematic errors in an experiment measuring four cross 
sections may be significant; this holds even more if experiments at different places are involved. 

In this work (and in [5]) the effect of next-to- leading order corrections (NLOC) has not 
been considered. A number of other cases suggests that their effect on the asymmetries will 
be less important than that on the cross sections; a partial understanding can be found in 
Ref. [18]. With NLOC, the effect on the cross sections of changing the scales is, in general, 
reduced. Whether (and how much) this affect will be reduced on the asymmetries is unclear. 
Unclear also is to what extend NLOC affect the combinations A (tt) and A (K) , which at LO 
isolate the subprocess — ^ qq. Anyway, the interest in reaction (1.1) as a possible probe 
of Ag makes imperative the determination of NLOC. 
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FIGURE CAPTIONS 



Fig. 1. Results when each of the final hadrons Hi, i — 1,2, is tt"*" or tt". 

(a) Differential cross sections Ada / d(l)idx for direct and resolved 7 contributions for 
Q = Qc — kiT + k2T- A, B and C refer to the parton distributions of Ref. [9]. 

(b) Asymmetries A = {Ada / d(f)idx) / {da/d(j)idx) and their variation with changing the 
scales in the range Qc/2 < Q < 2Qc. Strong lines correspond to the scale Q = Qc- The 
bands with forward and backward slanted hatches show this variation for sets A and C 
correspondingly. For set B the variation is not shown. 

Fig. 2. The same as in Fig. 1 when each of the final hadrons Hi is or K~ . 

Fig. 3. The predicted asymmetries at = 7.18 GeV together with the recent experimental 
result of Hermes collaboration [7] . 

Fig. 4. Asymmetries A = [Ada / dclyidxr) / {da / d(f)idxT) ■ 

(a) When each of the hadrons Hj is 7r+ or tt". 

(b) When each of Hi is K+ or K'. 

Fig. 5. Asymmetries for the combinations (6.1) and (6.2) for cross sections Ada / d(j)idx. 

(a) For the combination (6.1) of pions 

(b) For the combination (6.2) of kaons 

Fig. 6. The same as in Fig. 4, but referring to the cross sections Ada/d(f)idxT 
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